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Abstract 

We introduce a generalization of the extended Airy kernel with two 
sets of real parameters. We show that this kernel arises in the edge 
scaling limit of correlation kernels of determinantal processes related to 
a directed percolation model and to an ensemble of random matrices. 

1 Introduction and results 

The Airy kernel is one of the most fundamental objects of Random Matrix 
Theory. The determinantal random point process governed by the Airy ker- 
nel describes the behavior of the largest eigenvalues of large Gaussian Her- 
mitian random matrices (a.k.a. GUE - Gaussian Unitary Ensemble), see 
[6] , [TO] , [25] , last passage time in directed percolation models in a quadrant 
[T7] , asymptotics of the longest increasing subsequences of random permu- 
tations [I], and it also appears in many other problems whose list is too long 
to be included here. 

The Airy kernel has a time-dependent version usually referred to as 
the extended Airy kernel. Originally obtained in [23] via asymptotics of 
a polynuclear growth model in 1+1 dimensions, the extended Airy kernel 
arises in virtually every problem where the usual Airy kernel comes up, pro- 
vided that the probability measure in question is equipped with a natural 
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Markov dynamics that preserves the measure. In particular, it describes the 
edge scaling limit of Dyson's Brownian motion on GUE, the change in the 
quadrant last passage time when the observation point moves [18], and edge 
behavior of large random partitions under the Plancherel dynamics related 
to the longest increasing subsequences of random permutations [3J. The 
Extended Airy kernel has also appeared in a much earlier note [19], and we 
are very grateful to a referee for pointing this out. 

In [2] it was demonstrated that the Airy kernel is not stable in the sense 
that it can be naturally viewed as a point in a family of Airy-like kernels 
indexed by a finite set of real numbers. In the context of GUE, those are 
the (scaled) eigenvalues of a deterministic perturbation of finite rank, and in 
the percolation context the parameters correspond to a few defective rows 
or columns in the quadrant. For Wishart ensembles of random matrices, the 
same family of kernels was later obtained in [7]. Time-dependent extensions 
of kernels from this family appeared in the recent work [14] on asymptotics 
of the totally asymmetric simple exclusion process (TASEP). 

The main goal of this note is to introduce a new Airy-like time-dependent 
correlation kernel with two sets of real parameters. We obtain it as a limit 
of a directed percolation in a quadrant which has both defective rows and 
columns. We also show that the "static" version of the kernel arises in the 
edge scaling limit of a certain Wishart-like ensemble of random matrices. 
We were unable to obtain the extended version via random matrices but we 
do believe that it should be possible. Our kernel generalizes all the kernels 
mentioned above. 

Let us describe our results in more detail. 

Consider a directed percolation model with exponential waiting times 
defined as follows. Let tti, . . . , tt p , ttx, . . . ,tv p be fixed real numbers such that 
7Tj + 7Tj > for any 1 < i, j < p. Let W = (Wij)i,j=i,..., P be a p x p array 
of independent exponential random variables with E(Wi,-) = (7Tj + vrj) -1 . 
For any 1 < N < p, we consider the so-called last passage time in this 
percolation model: 

Y(N,p):= max £ W ij; (1) 
(ij)eP 

where II is the set of up-right paths from (1, 1) to (N,p). The random vari- 
able Y(N,p) has a natural interpretation in terms of queuing theory, due to 
the result of [13] . This is the exit time of the pth customer in a series of N 
files where the service times Wij depend on both the file and the customer. 
This random variable also has an interpretation in terms of TASEP, which 
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is a model of interacting particles on Z. One starts with the initial config- 
uration rjo(i) = meaning that only the negative sites are occupied. 
Then, if the site i + 1 is unoccupied, the particle at site i jumps to site i + 1 
after a random waiting time. The waiting times are independent exponential 
random variables whose parameters depend on the particle and the number 
of jumps already performed by this particle. One can think of Wij as of the 
waiting time of ith particle from the right performing the jump number j. 

We first prove the following result. 

Let Xtv be a p x N random matrix with independent complex Gaussian 
entries 



Theorem 1. Let Ai be the largest eigenvalue of X^X^. Then, for any x, 



A natural question is then to investigate the above connection as a pro- 
cess. Consider a sequence of growing random matrices (-Xjfc)jfc=i,...,p where 
Xk + i is obtained from AQ. by adding one column with random Gaussian 
entries (with the appropriate variance). We can then consider the joint dis- 
tribution of the largest eigenvalues of the random matrices X/.X?, 1 < k < p. 
Simultaneously, one can consider the joint distribution of the random vari- 
ables Y(k,p), 1 < k < p. Are these joint distributions the same? We cannot 
establish that the equality actually holds, due to the fact that the compu- 
tation of the joint eigenvalue distribution of the random matrices X/.X? is 
not an easy task0 Nevertheless, we can study a determinantal point process 
which occurs naturally in both models and obtain a new limiting correlation 
kernel, which generalizes the extended Airy kernel. 

Let Ji, J 2 be given integers, and X = {x 1 ,x 2 ,.. . ,x J:L }, Y = {y 1 ,y 2 , ■ ■ -,yj 2 } 
be given sets of real numbers satisfying Xj > y* for any 1 < i < J\ and any 
1 < j < J 2 . Let 7 and V be the contours defined on Figure Q] below. 
Denote by KAi{ti-,x;t 2 -,y) the extended Airy kernel 



1 As was pointed out to us by Peter Forrester, the equality can be established in the 
degenerate case when all Hi tend to the same constant using the techniques of [12], see 
Appendix to [TT] . 





¥(Y(N,p) < x) =P(Ai < x). 




(3) 
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Figure 1: The contours of the parameterized extended Airy kernel. 



Definition 1. The extended Airy kernel with two sets of parameters is de- 
fined by 

KAi;X,Y(h,x;t 2 ,y) = K A i(h, x; t 2 , y) + 7^~^2 J da J dr 



T + t 2 -a-ti V 11 ti + a- Xi Al t 2 + r-j/< / ' W 

\Z = 1 2 = 1 / 

where the two contours 7 and T are chosen as on Figure U\ 

Remark 1. For any fixed values of ii and t 2 , the kernel of Definition [T] 
is a finite rank perturbation of the extended Airy kernel. Indeed, by a 
straightforward computation, one finds that 

KAi-,x,Y{ti,x;t2,y) - K A i(ti,x;t 2 ,y) 

J1+J2 



f da f dTe y^ 3 /3-™+° 3 /3 



y y ^-i)^ {T + t2 _ a _ tl) k 1+ k 2 -i ^ 

l<h<-<i kl <Jll<h<-<j H <J2 nfcl(^ + *l -Xil)T\lLl(T + t 2 - yjll ) 

Remark 2. One does not have to stop at considering finitely many pertur- 
bation parameters. By taking limits with number of Xj's and y/s going to 
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infinity, one arrives at the following kernel. Let {af and {bf } c *i 1 be four 
sequences of nonnegative numbers such that X^i( a f + bf) < oo, and let 
be two positive numbers. Set 

$ . (z) - e c+z+c- z ~i TT (l+bfzKl+bTz- 1 ) 
aM{ ] ~ fi(l-a+,)(l-a-,-l)- (6) 

Then the kernel 

KAi;a,b,c(h,x;h,y) = K Ai(h,x;t 2 ,y) 

+ J l_U\ dT e^W^ fjatefelM-! 

where the contours are chosen so that all points 1/af — t\ and — t± are 
to the right of 7, and all points —1/bf — t% and — b~ — t 2 are to the left 
of r, is readily seen to be a limit of kernels of Definition [U Interestingly 
enough, functions © also parameterize stationary extensions of the discrete 
sine kernel, see [3]. They also appear as generating functions of totally 
positive doubly infinite sequences [8], and as indecomposable characters of 
the infinite-dimensional unitary group, see |22j and references therein. 

The main result of this note is the proof that the kernel KavX,y arises 
as a scaling limit of correlation kernels of the determinantal point processes 
related to the directed percolation model in a quadrant, defined as follows. 
Let < t < 1 be a given real number and assume that there exist two 
integers J\, J 2 , real numbers X{, 1 < i < J%, yj,l < j < J 2 independent of p 
such that: 

^ , Xi -1 7 1^7 

TT i = — H -pr, 1 = 1, ...,Jl\ TTi = l,l> Jl', 

y/t + l ap 1 ' 6 

Vi yj . . . 

TTj = -= J -pr, J = 1,...,J 2 , 7T j = 0,J>J 2 . 

We consider the associated directed percolation model ([1]) and prove the 

following result. 

Set 

,(t(i + Vi)) 2/3 



a s = t + 2 -7:; s, Y s = Y (pa s ,p) and a s = p (1 + y/a~ s ) . 

pi- 1 6 

Let also A e s xt be the stochastic process whose finite dimensional distributions 
are given by 

P(A t e f <&,..., < £m) = det(J - /K Ai;Xi y/) L2(rmxRm) , 
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where f(tj,x) = l x >^ and T m = {t\, . . . ,t m }. Here T m x K m is equipped 
with the product measure d[i (8> dX where d\i is the counting measure and dX 
denotes the Lebesgue measure. 

Theorem 2. As p — > oo, 

p-l/3^/6^ + ^t)4/3(y s _ ^ A ext^ 

in the sense of convergence of finite dimensional distributions. 

The paper is organized as follows. In Sections [2] and El we study the 
directed percolation model and its connection to the random matrix model 
defined in ((2|). This connection is established using the Robinson- Schensted- 
Knuth correspondence. In Section HJ we study the correlation function of 
the point process induced by the joint distribution of {Y(k,p)}i<k< p - Last, 
in Section we consider asymptotics of these correlation functions thus 
obtaining the kernel of Definition [TJ 

Acknowledgments. The first named author was partially supported by 
the NSF grants DMS-0402047 and DMS-0707163. The second named author 
was visiting UC Davis when this paper was written. She would like to thank 
especially A. Soshnikov who allowed her to spend a great year at UC Davis 
and encouraged this work. 

The authors are also very grateful to Peter Forrester for a number of 
valuable remarks. 

2 Last passage time in percolation models with 
exponential waiting times. 

We start with some reminders. We denote a partition by A through- 
out the paper. A partition is an infinite sequence of non-negative integers 
(Ai, A2, • • • , Ajv, . . .), where Ai > A2 > • • • > Xn > • • • , with finitely many 
nonzero entries. The Schur measure introduced in |21j is a measure on parti- 
tions which assigns to a given partition A a weight as follows. Let a\, . . . , ajy 
and &i, . . . , bp be given nonnegative real numbers. In what follows we assume 
that dibj < 1 for any 

Definition 2. The Schur measure with parameters a = (ai,...,a/v) an d 
b = (b\, . . . ,b p ) is a probability measure M on partitions assigning to a 
partition X the weight 

M(X) = L x (a)s x (b). (7) 
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Here Z is a constant, and s\ denotes the Schur symmetric functions param- 
eterized by A. 

The normalizing constant Z is computed by the well-known Cauchy iden- 
tity for Schur functions, which implies Z = j j (1 — ciibj). 

l<i<P,l<j<N 

The Schur measure naturally occurs in some directed percolation models, 
as we now recall. Let W = (iUy)i<i<p l<j<iV be a p x N random matrix 
with independent entries with the geometric distribution: 

F(wij = n) = (1 — aibj)(aibj) n , for any n = 0, 1, 2, . . . . 

One can associate to this random matrix the so-called last passage time in 
a directed percolation model defined by 

L{N,p) := max ^ w {j , (8) 

where II is the set of up-right paths from (0,0) to (p,N). This last passage 
time can also be understood as the exit time of the iVth customer in a series 
of p files with independent geometric waiting times of expectation depending 
on both the files and customer. 

The distribution of the random variable ([8]) can be conveniently ex- 
pressed in terms of the Schur measure, as observed by K. Johansson in 
[TT] . see also [18]. He showed that 

¥(L(N,p) <n) = ^Y, s a(«>a(&). (9) 

Ai<n 

This is a corollary of a more general fact that the Schur measure with pa- 
rameters a and b is the image of the random integer valued matrix W under 
the Robinson-Schensted-Knuth correspondence, see [17], [15] . 

We now turn to the description of a continuous version of the Schur mea- 
sure. Another treatment of the same object can be found in Appendix A to 
[12]. Set 

a,i = 1- j-, i = l,...,p, and bj = 1 - j = 1,...,N. 

It is quite clear that the distribution random variable Yl(N,p) := j^L(N,p) 
should converge as L — > oo to that of the last passage time in a percolation 
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model with i.i.d. exponential random variables with expectations 1 / (vr^ ) • 
Define Y(N,p) as in ([I]). Taking the limit L — » oo in ([9]) and denoting 
x i = Ai/L, one readily obtains in the case where N = p that 

i r p 
F(Y(p,p)<x) = — det(e-^)J det (e^Y jjdxi, (10) 

where Z PjP = det ^.^. ^ and I = [0, x]. The probability measure 
defined by the density 

' det(e- Wia, i)J. =1 det(e-* < ^) P (11) 



is the continuous version of the Schur measure. 

The case N < p can be handled via the limit transition bj — > 0, j 
N + 1, . . . ,p, from the case N = p. The analog of (fTUj) reads 

¥(Y(N,p)<x) = J- / det(/ i (x i ))f detfe-^Y TTdx,, (12) 
^V.p J IN > J x 1 1,0=1 j-J; 

with 

1 f u k ~ 1 e~ ux du 



and the integration contour going around the poles 7T±, . . . , tt p . 



The expression det (fi(xj))^- =l can be obtained as the limit of the Jacobi- 



Trudi formula for s\(a); it is also the limit, up to a constant, of the ratio 

)i,j=l 



det (e-^Oi 



11 Oi - xj) 

N+l<i<j<p 

as xn+i, ■ ■ ■ ,x p converge to 0. 

3 A random matrix model associated to the con- 
tinuous version of the Schur measure. 

The goal of this section is to prove Theorem [TJ 

Let Xn be a p x N random matrix as in ([2]). Set then Mn = X^X^. 
Thus defined random matrix ensemble is a natural generalization of the 
much studied complex Wishart ensemble. In what follows we call it the 
generalized Wishart ensemble. We show that the probability distribution of 
the largest eigenvalue of has a density given by the right-hand side of 

(nrj. 
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Let us first consider the case where p = N . Similarly to the case of the 
ordinary Wishart ensemble, one sees that the generalized Wishart ensemble 
is defined by the probability density with respect to Lebesgue measure L p 
on the space of complex matrices of size p x p: 

dF j^ = constp exp {-TrSiX p X* p - TrS 2 X*X p }, (13) 
(Llj p 

where S\ = diag(7Ti, . . . , n p ), S 2 = diag(7ri, . . . , tt p ), and const p is a positive 
constant. Just as for complex Wishart ensembles with non identity covari- 
ance matrix, the joint eigenvalue distribution of the generalized Wishart 
ensembles can be explicitly computed. 

Denote by x\ > x 2 > ■ ■ ■ > x p the ordered eigenvalues of the sample 
covariance matrix X p X*. Note that these eigenvalues are also the squared 
singular values of X p . Let f(xi, . . . , x p ) denote the density with respect 
to Lebesgue measure of the joint eigenvalue distribution induced by the 
generalized Wishart ensemble. 

Proposition 1. One has 

f(xi, . . . ,x p ) = — det (e-^Y . , det (e-*i x "Y . 

Z P)P K h >3 =1 V Jj,k=l 

Proof of Proposition Q3 Introduce the polar decomposition of the p x p 
matrix X p : One has 

X = UDV with U G U(p), D = diag(^/il, . . . , ^), and V G U(p). 



The joint eigenvalue distribution induced by the probability measure f)13j) 
can now be computed thanks to the celebrated Itzykson-Zuber-Harisch- 
Chandra (IZHC) integral. We have 

/(*!,... ,x p ) = const -V{xf [ e- T ^ UD2u *dU [ e~^ VD2y ' dV. 

JV{p) JV{p) 

In the above expression, V(x) is the Vandermonde determinant: V{x) = 
Y\i<ji x -i ~ x j)- The IZHC formula yields 



/ 



p -Tr5it/D 2 C7* r n _ det ( e " ' " >i,j=l n 

The proof that the largest eigenvalue of X^X^ has the same distribution 
as the random variable Y(N,p) can now be obtained from Proposition Q] and 
formulas of the previous section by the limit transition 7Tjv+i, . . . , n p — > 00. 

This finishes the proof of Theorem [TJ □ 
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4 A continuous version of the Schur Process 



In this section, we define a random point process which is a continuous 
version of the Schur process introduced in [23] . We consider the probability 
distribution on nf=i ^+ w hh density w.r.t. Lebesgue measure given by 

^det (e-^) P fldet ( e -*** + i(*? +1 -*< fc >) e ~-< (14) 

Here we used the convention that x k+1 = for any 1 < k < p — 1, and the 
notation 



e 



* Hx - y) = le- Hx ~ V \ *>y, 

I 0, otherwise. 



The probability distribution (|14p naturally arises here since the distribution 
of maXj = i j ... )P {z?} is equal to the probability density function of Y(p,p) (see 
Formula ([IU])). 

Let C (resp. C') be a contour encircling the {^j}j=i,...,p (resp. {— tTj}^!,...^) 
such that the two contours do not cross or contain each other. Set 

* r , s (u,v) = l r<s l u<v / e^-") TT ^—dw. 



The main result of this section is the following statement. 

Theorem 3. The random point process on {1, . . . ,p} x R defined by the 
density [14\ ) is determinants, and its correlation kernel has the form 



K(r, u; s, v) 

= TT-y2 ( f dz f dw tv ) 7^ 11 *r,.(«,«). 

(2m) Z Jc J C , W - Z [[i =1 {W + TTl) Z - 7Ti 

(15) 

Let us briefly discuss the connection with the (discrete) Schur process. A 
version of the Schur process has a natural interpretation in terms of the last 
passage percolation model discussed in Section [2j Let again W = W(p) be 
the p x p matrix filled with geometrically distributed integers, and let W^ k \ 
k < p, be the px k matrix made of first k columns of W. Denote by A( fc ) the 
image of under the Robinson-Schensted-Knuth correspondence. Then 
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the joint distribution of (A' 1 ', . . . A^) is given by the Schur process: It has 
the form 

const ■ s X ( P )(ai, . . . , a p ) s X ( p) / X ( P -i) (6 p )s A ( P -i)/ A ( P -2) (bp-i) • • • s A (2)/ A (i) (&i), 

(16) 

where the notation s a a, stands for the skew Schur function. Considering 
the case where the Oj's and bj's approach 1 (m = 1 — vrj/L, foj = 1 — 
ttj/L and L —* oo), we can then define the continuous limit of the Schur 
process, which leads to (|14p . In the context of the percolation model, the 
probability distribution (JT1J) can be understood as the joint distribution of 
the random Young diagrams obtained by the RSK algorithm applied to 
matrices filled with independent but not identically distributed exponential 
random variables; the expectation of the (i, j)th entry is equal to (7Tj+7rj) _1 . 

Theorem [3] could be derived from a limiting argument for the correlation 
kernel of the Schur process, but we prefer to give a self-contained random 
matrix oriented proof of Theorem [3] below. 

The probability distribution (|14ft can also be viewed in the random ma- 
trix theory context of the previous section. For 1 < k < p, define X^ to 
be the p x k matrix whose k columns are the first k columns of X p . Then, 
Mfc = Xf,Xt is a p x p random matrix of rank k. We denote by x*, 1 < % < k 
its nonzero eigenvalues. The formula (|14p provides a good candidate for the 
joint distribution density of {x^} in the sense that its projections to {xf^} 
with fixed k coincide with the densities of eigenvalues of M^. Although we 
were unable to verify that (I14h is indeed the joint eigenvalue density for 
(Mi, . . . , M p ), see the footnote on page 3 for a reference to a partial result. 

Proof of Theorem [3]. We first consider an (algebraically simpler) aux- 
iliary distribution and then use an appropriate limit transition to compute 
the correlation functions associated to (fT4"|) . 



Instead of ([T4"j) let us consider the probability distribution defined as fol- 
lows. Let < Ti < T2 < • • • < Tp_i < T p be positive numbers, and T% = n±. 
Define a probability distribution on (M^j_) p = {x^}k,i=i...,p by the density 




(17) 
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where fa^x^x) = e * T * X , (f) r>r+ i(x,y) = e **r+i(w x) and ^ pp+1 (x, xp +1 ) 
e^* Set 



^oA x h v )= / 00,1(4^1) 

's-2 

Y\ 4>k,k+i{x k ,x k+1 ) J s _i )S (x s -i, u) \\ dxi 

Jfc=l 

*r,p+l («, xl+x ) = / <^r,r+l(«,^r+l) 
jR p + - r 

/ P-1 \ P 

n ^.fc+i^Aj^+i) ^p.p+i^p, ^+1) n dxi, 



s-l 

s-2 \ s-1 

vfc=l / i=l 



\fc=r'+l / i=r+l 

^r iS (u,v) = l r <s 4>r,r+l(u, X r ) 

(s-2 \ s-l 

]^[ 0fc,fc+i(a?jfc,£fc+i) s _i )S (a; s _i,i;) JJd^. 
ft=r+l / i=r 

Lemma 1. For any 1 < i,j < p and u, u > we /zaue 

^r„(xn,f) = & — TT — dw, 

°' sV °' ; 2mJw + T l £} 2 w + fc k 

p 



^ )P+1 («,4 +1 ) = e-^ f[ -i 



Aij :- 



TTj + 7T fc 

/ ^o,s( x o^ u )^r, P +i(u,^j)du = — — — TT — — -. 
*r,.(«,«) = lr<»l«<«^T /e^-") TT — 

k=r+l 

The integration contours are positively oriented loops that contain all poles 
of the integrands. 

The proof of Lemma Q] consists of inductions on s for s (xq, v), on 
(jp—r) for ^ riP+ i(n, x J p+1 ), and on (s — r) for \& rjS (n, u). The formula for 
is proved by a straightforward residue computation. 

Let us now apply the Eynard-Mehta theorem (see [§] , [20J , [26] , [IB] , [5] ) 
to compute the correlation functions of (|17p . For 1 < r, s < p, denote 

^J,s( U > V ) = l r<s I 4> r ,r+l(u,X r ) 



k=r+l 

-oo 
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(s-2 \ s-1 

Yl 4>k,k+l{xk,Xk+l) J 4>s-l,s{x s -i,v)Y\_dXi. 
k=r+l / i=r 

Proposition 2. The random point process on {1, . . . ,p} x M + defined by 
the measure |i7| ] is determinantal, and its correlation kernel can be written 
in the form 



Kx(r,u; s,v) = — — d> dz d> dw 



{2iri) 2 J Cl 7 C , w-z 



LU=2(™ + **) f=\ ( z - n)(vi + Ti) 

where the contour C\ encircles the TTj,j = 1, • • • ,p, the contour C[ encircles 
the —TTj, —Tj, for j = 1, . . . ,p, and the two contours do not cross or contain 
each other. 



Proof of Proposition [2} The Eynard-Mehta theorem implies that the 
random point process in question is determinantal, and that its correlation 
kernel can be written as 

p 

K T (r,u;s,v) = ^ ^r.p+i^.^+^^Wi^-^vK")- 

Using the formula for the determinant of the Cauchy matrix one explicitly 
computes A^- , which together with the formula for from Lemma [T] 
yields 

K T (r, u; s, v) + ^f r ,s(u, v) = 

e*- + i(«,4 + i)^ / e- vup : ?\ n ^ n — *». 

^ p+ 2vrz / Uk=2( w + Tfc) f = i w + T 3 y t< - 

where the contour contains all poles of the integrand. A final residue com- 
putation yields the integral expression (|18p . □ 



We can now come back to the computation of correlation functions for 
the probability distribution (I14j) and finish the proof of Theorem El The 
probability distribution (|14p can be obtained from (|17|) by taking the limit 
T p > T p _i > • • • > T2 — > 00. The proof of Theorem [3] is now a straightfor- 
ward corollary of Proposition [2l □ 
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5 An extension of the Airy point process 



In this section, we first consider the case where TTi = and 7Tj = 1 for any 
i = 1, . . . ,p. We then show that in this case, the suitably rescaled correlation 
functions converge in the large-p-limit to those defined by the extended Airy 
kernel. Then, to define a new extended Airy-type kernel with parameters, 
we will allow a certain number of these parameters to depend on p and study 
the rescaled correlation functions. 



5.1 The simple case 7Tj = and 7Tj = 1 for any i: the extended 
Airy kernel. 

The extended Airy kernel is an extension of the well-known Airy kernel; 
it occurs for example as the limiting correlation kernel for the process of 
largest eigenvalues of Dyson's Brownian Motion on Hermitian matrices. 

Lemmas 2 and 3 proved in this section are also a part of Proposition 5 
of [IT]. 

Let KAi(ti, x;t2,y) be the extended Airy kernel defined in ([3|). Here we 
show that, when suitably rescaled, the asymptotics of the correlation kernel 

gpwv— pzu 



K p (r,u;s,v) 



P 



dz (b dw- 



w — z 



(27Ti) 2 J Cl J C , 



which is a rescaled version of (|15p . is given by the extended Airy kernel. 
Due to the choice of the n^s, and 7Tj's, one can write that 

K p (r,u;s,v) : = K p (r,u; s,v) + p^ r)S (pu,pv) 



(2ivr) s 



dz 



dw e P {F v>s (w)-F u , r (z)) 1 

w — z 



where F u>r {z) = uz + ln(z — 1) — ^ Inz. 

Let < t < 1 be some given real number independent of p. Let also ti, £2 
be given. In the following, we set 



tp + p 2/3 -^~ ^ Z-t-tl 



a 



tp + p 



a 



2/3 2Vt(l + Vt) 2 



a 



t-2 



SlP, 
S2P, 



£V6 



2/3' 



(l + ^) 2 + 



ay 

^2/3 



(20) 
(21) 
(22) 
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Here [x] stands for the integral part of x £ M. 

We first consider the case where s < r. Set V := {te^ 2m ^, t G K+} to 
be a contour oriented from bottom to top and 7 := {te^ tn ^ 3 , t € M + } to be 
oriented from top to bottom. 

Lemma 2. With the above rescaling, for s < r and zq := we have 

lim " e p F*A*)-FvA*>))k ( r , u- s, v) 

p—>OG p*/ «J 

e yt 2 — a;ti + § (if — 1|) /• gyi'-^ — rrs'+^- 



(2vri) 2 J J T t'-s' + t 2 -t 1 



/»oo 

/ e - A( * 1 "* 2) ^(y + A)yli(a; + A)dA. 
•/ 



Proof of Lemma [21 It is convenient to define u Q = (1 + \/^i) 2 an d 
(1 + y^) 2 - The reason for the above rescaling ([2"0j) - ([2"2"|) is that 

FuA z ) = F Uo ,r{ z )( 1 + °( 1 )) 

= (u z + ln(* - 1) - 8! lnz)(l + o(l)) := / ai (*)(l + o(l)), 

where the function f si (z) has a degenerate critical point at 

zo + Oip- 1 / 3 ). 



1 + 

In particular, one has that: 

-2(1 + ^1) 



f' sl ( z c)=fsM=0, and/-(z c ) 
To obtain the leading term in the asymptotic expansion of K p , we define 



the following contours. Using the notation w c = j^—^=, set 



C M = {z c + s'e ± ^/ 3 , < s' < 5 }, C[ A = {w c + t'e ±i27T ^, < t' < tfj}, 

where 5 D > and <5i > are constants that will be determined in the sequel. 
These contours are completed as follows. Set O = arg(z c + e) e i7r / 3 ) (resp. 
Q\ = arg(w c + 5ie 2tn / 3 )), where arg denotes the argument of a complex 
number. Set 

C[ 2 = {\w c + 5ie 2m/3 |e ie , 9 1 < 9 < 2vr - 61]}, 
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Figure 2: The two contours C\ and 



Ci,2 = {1 + \z c + <W /3 - l|e ie , \e\ < e Q }. 

The constant S Q is also large enough so that the z-contour encircles all the 
7Tj's (even in the case where some of them differ from 1). It can then be 
checked that 

±» ( r r«, + \\ - s2 ( 1 + v / ^)( s2 - s ( 1 -v / ^) + 2 v / ^) < n 

5s I 1 + \/*2 / " 2( S 2 + S + 1)( S 2_ S ^ +S2 ) 



for any s > provided > 5 — V96/2. If y/~S2 < 5 — V96/2, we set Si to 
be the smallest positive root of the polynomial X 2 — X(l — y/s2) + 2^/s2■ 
Otherwise Si is arbitrarily large. Similarly 

±»(f (? 4. seW3 ^ _ ^(i + v / ^)(^ 2 + Ki-v / ^) + 2 v ^i) 

for any s > 0. We also have that if \S\ is bounded 



Pfs 2 (w c + - pf S2 (w c ) - f£\w c )^ 



< supnl/ffHI < ^ <5 4 



Here u; = tt; c + lies in a compact subset $7ofC\{0,l}. 
To complete the proof one needs to verify that on remaining parts of the 
contours the integrand becomes exponentially small as p gets large. If we 
set w := \ w c + 5ie 2 " r / 3 |e* e , and using the fact that v Q = \w c — 1\~ 2 , one can 
check that 

* (^/« (")) = ~^ w) <S - k^Tp) < °' 
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if 9 £ [9 1 , tt] . It is an easy computation to check the remaining parts of the 
contours, and we omit it. 

If we assume that ti,t2,x,y lie in a compact set of R, then the above 
estimates imply that 

• It is enough to integrate over a neighborhood of radius p~ l / 3 p l / l2 ~ e of 
the critical points. 

• Inside such a neighborhood the Taylor expansion holds. 

Writing out this expansion explicitly yields Lemma [2j □ 
In the case s > r we cannot make the integration contours go through 
their corresponding critical points so that they do not intersect. We then 
modify the contours in a neighborhood of width p -1 / 3 of z c and w c so that 
the w— contour remains on the left of the z— contour. This does not modify 
the saddle point argument. We just need to consider the function ^f r ^ s 
separately, which is done below. 

Due to the rescaling of the correlation kernel, in the case where r < s, we 
need to consider the asymptotics of 



a ^(Fu.Az^-Fv.Aza)) P I c pw(v-u) dw 
p 2 / 3 2m J lo w s ~ r ' 

where "y Q is a contour encircling the pole w = 0. 

Lemma 3. For s > r, using the scaling \20\ ) -l23 1) one has 

li m " f ?F«A*o)-pFv,.(*o)K t , 

p— >oo p^/d 



ft?- th ro 



-e a " s / e 



-K*y-ti)Ai(y + X)Ai(x + X)d\. 



Proof of Lemma [3} Setting [3 = 2\/i(l + Vt) 2 /a, one has 



a c {pF u . r (*a)-pF„.M\ P I c pw(v-u) dw 
p 2 / 3 2ni jL o w s ~ r 



n 



,{pF v , r (zo)-pF v , a (z )} P 



dw- 1 ^ J 



p 2 /3 2TTi J 7 w /3p 2 / 3 (t2-ti)+0(l) 

4t 3/2 



x exp { p 1 / 3 I --^{t\ - t\)w + a{y - x)w + o(l) 
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Consider 

F{w) = (3{t 2 - h) 1+ /* w - (3(t 2 - h) Inw. 

It is not hard to see that the critical point of this function is wo = zq = 

and F"(zq) = (3{t 2 — ti)zQ 2 . A contour which satisfies the saddle point 
analysis requirement can be chosen as follows: 7o = 7i U 71 where 71 = 
{z + it, \t\ < z } U {zoe^+O^O < 6 < 3vr/4}. 
We obtain 



Oi P f 

lim -—j-exp{pF ur (z )-pF vs (z )}— (t 

p-*COp z /3 2iri J-y 



c pw(v-u) dw 



111 

'70 w 

a 1 /• rn»-^ — 4^72— ) s' 2 



1 f (s'(a(y-x) i( 
lim : — ^=^= / exp < — ^=^= — H \ds 



P-00 27Ti y/F»fa) Jm 1 V^o) 2 

= ^^= / -pf*", fa-^'-^-^'l.,, 

2^ v / 2(i 2 - ti) Jm \ 2 ^2(t 2 -ti) J 

1 f (^-x + tf-tl) 2 
= = exp < 

V^Mi) I - *i 

Proposition 2.3 of [IB] completes the proof. □ 



5.2 Extended Airy kernel with two sets of parameters 

We now consider the case where some of the n^s (resp. TVj's) differ from 1 
(resp. 0). This allows us to obtain a new extended Airy type kernel with 
two sets of parameters and prove Theorem 2. 

We assume that ([7]) holds true and that all the Xj's and y^s lie in a fixed 
compact set of R. We also assume that X{ — yj > for any so that the 
joint distribution (|14p is well defined. 

Theorem 4. With the above rescaling, one has 

lim -^-e pK > r{zo) - pF "^ Zo) KJr,u;s,v) = Ol) 

p— »oo p'/ 1 — ' 

where the integration contours 7 and T are chosen as in Figure [7J 
Remark 3. Theorem [3] readily implies Theorem [2j 
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Proof of Theorem [4} The proof relies on the same saddle point analysis 
of the correlation kernel as in the previous section. In the expression (|19p 
we replace 

W ~ Z Uk=l( W + *k) Z ~ T< 

by 

rifc=j 2 +i( z +^) -rj w - I i (rifcii ip+^t n»=i f=^f ~ i ) \ 

Tlfc= J 2 + l (™ + *k) ^-VTi^-Z J 

and observe that the second summand has no singularity at z = w. This 
allows us to use the same contour deformation as in the previous section, 
which directly leads to the result. □ 
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